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Abstract 

In this paper we will present a quantization method for SPP (Surface Plasmon Polariton) based 
on Green’s tensor method, which is applied usually for quantization of EM field in various dielectric 
media. This method will be applied for a semi-infinite structure, which contains metal and dielectric 
regions with one interface. Moreover, by introducing the quantized SPP, we will investigate the 
SPP propagation in the attenuating and amplifying systems. We will also consider two modes of 
SPP, i.e., coherent and squeezed states, and finally compare the propagation of these modes in the 
amplifying media. 

PACS: 42.50.-p 
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I. INTRODUCTION 


A SPP is a an electromagnetic 2-dimensional excitation existing on the metal-dielectric 
interface, which decays exponentially with diffusing depth in the metal. The term ’’surface 
plasmon” refers to the collective electrons oscillation on the metal surface , and ’’polariton”, 
which combine exciton and photon properties as a quasi-particle, indicates that the SPP 
involves coupling of electromagnetic waves and the dielectric excitations. Due to the unique 
properties of the SPP, such as capacity of confining the electromagnetic held [1, 2], it has 
become one of the interesting held of applied research [3]. 

By discovering the quantum behavior of SPP in theoretical and experimental investiga¬ 
tions [4-7], the possibility of using them in quantum technology, such as Quantum Infor¬ 
mation Processing (QIP) [8-10] and generating and controlling plasmonic excitation at the 
quantum level, has been provided [11, 12], 

In order to describe the plasmonic quantum mechanically, a quantization method of SPP 
is introduced in [13]. This formalism is based on Hopheld approach and therefore do not 
consider the dissipation within the medium. It is well known that the real part of permittivity 
is related, by the Kramers-Kronig relation, to the stored energy or dispersion within medium, 
and imaginary part yields the dissipation (or loss) of the energy within it. For lossy and 
dispersive dielectric, Huttner and Barnett have proposed a microscopic quantization method 
[14] . In their method the absorptive medium is assumed to be connected to a reservoir and 
represented by a collection of interacting matter fields. 

As another most interesting scheme for quantization, the Greens function method have 
been proposed and used by many authors [15-18]. The Greens function is constructed for 
classical Maxwell equations and contains, as a whole, the system and the source of the noise 
associated with the absorption of radiation in it [19-22], 

In this contribution we will try to apply this method for quantization of SPP held. More¬ 
over, we apply this approach to propagation of the SPP in the amplifying media. Nowadays, 
because of the significant applications of this held, many efforts have been devoted to the 
SPP amplification [23-28]. The quantum approach not only enables us to take into account 
the SPP amplification, but also the behavior of SPP’s different states can be distinguished 

[29]- 

The paper is organized as follows. In section 2 we review the quantization of EM-held in 
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dielectric media. Regarding to the SPP propagation along the interface of metal-dielectric, 
we study a single flat interface system in section 3. In section 4 we apply this approach for 
investigating the SPP propagation in attenuating and amplifying system for coherent and 
squeezed states of SPP. A conclusions is given in sec 5. 


II. BASIC CONCEPTS 


In order to obtain a suitable form of the Maxwell equations for quantization procedure, 
we review the basic concepts in classical and quantum electrodynamics [15-18]. 

At first, it is useful to express electromagnetic fields in frequency domain, and separate them 
in negative and positive frequency components. For instance the electric field is written in 
the form, 

E(r, oj) = E + (r, oj) + E~(r, oj), ( 1 ) 


here 


1 f + °° 

E + (r, t) = / duE + (r, oj) exp(— iut) (2) 

v2vr Jo 

and the other vector fields are expressed analogously. Then, in a medium without free 
charges and currents, the classical maxwell equations for positive part in the frequency 
domain are given by, 


V ■ D + [r, oj) = 0, 

V ■ B + (r, oj) = 0, 

V x E + (r, oj) = iojB + {r, oj), 

V x H + (r, to) = —iuD + (r,co), 


where D + (r, uj) and B + (r, co) satisfy the constitutive relations: 

D + (r,uj) = e 0 e(r,u)E + (r,uj), 
B + (r,u>) = fiQ H + (r,u). 


( 3 ) 

( 4 ) 

( 5 ) 


where e(r,oj) is complex dielectric function, its real and imaginary parts are related by 
Kramers- Kronig relation. Also, e(r, uj) is related to complex refractive index, n(r, u), 

e(r,u) = [ n ( r , a ;)] 2 , ( 6 ) 


3 




( 7 ) 


e(r, -u) = e(r,u)*, 


n(r, —co) = n(r, co)*. 


( 8 ) 


As a general property for e(r,oo) it is, in the whole upper half of the complex frequency 
plane, an analytic function without zeros. So that, with respect to the causality, e(r, co) has 
asymptotic behavior for high frequency [30, 31], 


e(r, co) —* 1, for co —* oo 


(9) 


In quantizing the EM-vector fields and constitutive relations (4) and (5), we have to take 
into account the fluctuation- dissipation theorem. According to this theorem, an unavoid¬ 
able consequence of the radiation absorption (or dissipation) in reservoir, is injection the 
additional noise (like thermal radiation) from environment to the fields. The noise term can 
be expressed by introducing a polarization operator in Eq. (4), 


D + (r,v) = e 0 e(r,u)E + (r,ui) + Pj(r,w). 


( 10 ) 


So the quantized Maxwell equations are obtained as: 


V • [e 0 e(r,u)E + (r,uj)} = ptr(r,u), 

V ■ B + (r, uj) = 0, 

V x E + (r, oo) = iu)B + (r, co), 

V x 5 + (r,(n) = - iup 0 e 0 e(r,uj)E + (r,uj ) + p 0 j^(r,uj), 


( 11 ) 

( 12 ) 

(13) 

(14) 


here p%(r, u>) and j^(r,uj) are noise charge density and noise current density respectively 
which are related to (r, co) by, 


Piv(A^) = -V ■ 

Av(A^) = -iujP£(r,u). 


(15) 


It is suitable to obtain a partial differential equation for A(r, co). By expressing the held 
operators in terms of vector potential operator A(r, co): 


E + (r,co) = icoA + (r,cu), 
B + (r,co) = V x A + (r, co). 


( 16 ) 
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(17) 


and applying Eqs. (13) and (14), one obtains: 

to 2 

- V x V x A + {r,u:) + — e(r,oj)A + (r,uj) = ~n 0 j£(r,u). 

A standard method for solving this equation is using the Green’s tensor, 

/ +oo 

d r'G(r,r',uj).j^(r',uj), (18) 

-OO 

and therefore it has to solve the following equation, 

, u 2 

— V x V x G(r, r\oo) H—-e(r, oo)G(r, r , oj) = I5(r — r). (19) 

c 2 

where I is the unit tensor. The Green’s tensor has some general properties, which are given 
in appendix A. In the held quantization one requires the vector fields satisfying the canonical 
commutation relations at any point of space, 

[A( r , t), —e 0 Ej(r,t)] = iH6(r-r')6ij. (20) 

The main reason for accuracy of this relation in absorbing media is the noise sources. The 
commutation relation for the noise current operators is given by: 

[ilv +) (G w )Ow VV)] = HAI 6(r-r)6(u}-u), 

= 0 . ( 21 ) 

where a(co) is determined such that Eqs. (20) and (21) be consistent with each other. On 
the other hand the noise current is related to Langevin operator, 

jN ] (r,uj) = sja(u)f(r,u ). 

Also one can obtain the propagators of noise current [15-17], 

(0 = (0| (r\uj')\0) = 0 

(0|j^“ ) (r,o;)j^ ) (r',a; , )|0) = 0 

(0|j^ ) (r,o;)jj7 ) (r',a; , )|0) = \a(u)\5(r - r')5(u - u) (22) 

III. QUANTIZATION OF THE SPP FIELD ON A METAL- DIELECTRIC IN¬ 
TERFACE 

In this section we consider the quantized held of SPP, which propagates along the metal- 
dielectric interface lying in xy-plane, as shown schematically in Fig.l. For simplicity we 
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FIG. 1. Semi-infinite media, SPP propagates along the interface of metal- dielectric. 

assume that the both media are homogeneous, so that the related dielectric constants are 
independent of space points. The dielectric constant corresponding to the whole space is 
written as, 

e(r, u) = e m (uj)Q(-z ) + e d (co)Q(z), (23) 

here e m (u>) and e d (u>) are dielectric functions of metal and dielectric media, respectively 
and 0 is the step function. In order to use the Eqs.(17) and (18) as starting points for 
quantization process, we must derive the Green’s tensor. 


A. Construction of the Green’s tensor 


For construction of the Green’s tensor we use the method of eigenmode expansion [32- 
34], In this approach, we consider the generalized form of Eq.(17), which contains a set of 
eigenmodes and eigenvalues (A n ,X n ), 


u 


V x V x A n (r) + ^e(r, w).4„(r) = e(r,w)A„Ai(r). 


(24) 


These eigenmodes satisfy the orthogonality condition: 


"- 1-00 


t(r,u)A„(r) • [AnWl'dV = 


(25) 


By obtaining the explicit form of eigenmodes and calculating the eigenvalues and normal¬ 
ization factor from Eqs. (24) and (25), the Green’s tensor is given by, 


G(r y) = ^ MdMT. 

An 


( 26 ) 
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By this method the Green’s tensor, for geometry shown in Fig 1 can be derived. The SPP’s 
eigenmodes can be written, 


(x, z) = (ix + ^z)e^e ik * x , z < 0 


A kx (x , z) = (ix - —z)e { ~ voz) e ik * x , z > 0 
G) 

here /y and n 0 are out of plane decay constant and are given by: 


(27) 


Vr=d—\k x \, 

Cd 


G> = \l —\k x 


(28) 


k x is the in plane propagation wave vector. The details of this calculation have been given 
in appendix A. By inserting Eq.(27) in to Eqs. (24) and (25), A n and N n are obtained: 


\ - b 2 b 2 €d + e ' n 

A n — fi-o 

^m^d 


Nn \k x \ 


~e m (e 2 d - 4J 


(29) 


According to Eqs. (26) and (27), for z, z < 0 and z > Q,z < 0, the Green’s tensor can be 
written as: 


G(r,r',u 0 ) = -C d k x { 


I k x 


(Mpp - k l - *o+) 


[(ix + — z)e^ UiZ ^Q(—z) + (ix — — z)d- I/ °^@(^)]x 
^ v 0 

[(ix + —z)e ( ' UiZ ^Q(-z') + (ix — —z)e l '~ U0Z )0(z)] 
G G) 


v ik x (x—x ) 


}. 


(30) 


here 0 is the step function and i0 + in the denominator is necessary for obtaining the retarded 
Green’s tensor and: 


C 


k 2 

^spp 


C-m^d 

7 .; w 

\Cd CmjTT\ / 

V Cd 

k 2 _Jn&_ = ( k ' + k y m 

£d + Gn 


(31) 

(32) 
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By using that lim e _> 0 - = PV —b iTrS(y) and the residu theorem one can evaluate the 

y-ie y 

integral in Eq.(30). This result in the far- held approximation yields: 

G(r,r',uio) = — iC x 

[(ix + ^z)e^ iz) Q(-z) + (ix-^z)e^ v ° z) Q(z))x 

Vi V 0 

[(ix + -^-z)e^ iZ )©(— z) + (ix — — -z)e^~ U0Z ) ©(^)] 

© v 0 

e ik spp \x-x\ ( 23 ) 


where C = ttC. 

B. Field quantization and Canonical commutation relations 

By substitution the Green’s tensor Eq. (33) into the Eq. (18), the vector potential 
operator is now obtained, 

i+(r , oj) = iy 0 C x 

[(ix + -^-z)e ( ' UiZ ' l Q(—z) + (ix — ^^-z)e^~ I/0Z ^Q(z)]x 
Vi v 0 

/»Od /»oo 


[dz dx e ik ^ x ~ x l j+(x,z',uj)- 


' —OO J — OO 


[(ix + ^z)e^G(-z') + (ix - k ^z)e^ oz) e(z)}. 

Vi Vq 


(34) 


The noise current contains two components, which are related to two regions (metal and 
dielectric), 


3n \ x , z, u) = j { N )m (x, z, u)Q(-z) + ] { N )d ( x , z, u)Q(z) 

= (sja m (oj)®(-z ) + y/c^(u})Q(z))f(x,z,uj) 


(35) 


It is convenient to separate the vector potential operator into two components, rightwards 
and leftwards, based on the propagation direction, 

ii+(r ’“ ) = i( ar ) ^ °c' x 


{(a + hs£i) e ("<’>e (- z ) + (it - *se*u-w>o 


ix - —z)e K - u " ) ©(^)}x 

vo 


[a R (x,uj) + a L (x,u)\, 


(36) 
















where Or(cil) is annihilation operator for rightwards (leftwards) SPP modes and have the 
explicit form, 


where 


a R {x,uj) = -^-2 
k 


dx dz e ikspp ( x ~ x )x 


r — OO «/ —OO 


[(z£ H——i)e^ iZ • j^ m (x, z, u)Q(— z)0(— 2 /)+ 

(ix — ^Lz)e^~ U0Z ^ ■ j^ d (x, z,u})Q(z)Q(z')], 

^0 


2fc” 


a L (x,o;) = (www ) 5 

(3(uj) 
k 


/»oo /»oo 

1 / / , / , / 


da/d^e-^^-^x 


[(z£ + -^ 2 )e (i/iZ ) • 2 ,a;)©(- 2 )@(— 

Vi. 


" o; ) ■ ^,w)@(s)0(s )], 

/3M = (i + + (1 + IM!)hM 

Kl 0 + ^ ror ^0 + ^0 


(zx-—i)e 

z'o 


(37) 


By this definition, one can obtain the bosonic commutation relation at the same positions, 


[a R (x,uj),a f R (x,uj')\ = [a L {x,u),a' L (x,iJ)\ =5{u-J), 
[a R {x,u),a\{x,J)\ = [a L (x,u),a} R (x,u)\ =0. 


(38) 


The representation of the vector potential in quantized scheme Eq.(36) is acceptable when 
the canonical commutation relation (20) is satisfied simultaneously. In order to show the 
accuracy of the relation (20), it is advantageous to derive the commutation relations for the 
creation and annihilation operators at any two points of space along the x direction, 


[a R (x,u),a ] R (x ,u)\ = 
5{oj — l 0 ) exp(zfc (x — 
[d R (x,u),d[(x ,u)] = 

S(u> — u )0(x — x 

k 


-- [d L (x',u'),d[(x,u)} = 

x )) exp {—k \x — x |), 

= [a L (x ,u),a' R (x,u)] = 

?? 

- exp (—k (x — x )) sin k (x 



(39) 


(40) 


The relations (39) and (40) are the same such as evaluated in Ref. [15] for field operator 
in dielectric medium. Matloob et al. in [15], have interpreted these relations, and also 
mentioned that the presence of the loss leads to coupling the rightwards and leftwards 
operators. By applying these relations, the canonical commutation relation is written in the 
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following form. 


[A(r,t),-e 0 E(r',t)}= I duj\C'\ 2 nl x 


I k 


spp 


JO 

|2 p ik app \x-x | e ~ik* pp \x-x' \ 


2k" n 


k 




k 


+ 


spp 


k 


-}x 


spp 


[(ix + n) + (fa: — ^^z)e ( ^©(z)] x 

[(—ix H— s -^-z)e^ z ®{—z) + (—fx —-^z)e _l '° z 0(/)]. 


ZA - 




(41) 


By some calculations, done in appendix B, the convenient form of the above equation can 
be achieved, 


[A(r,t),-e 0 E(r' ,t)\ 


do; 


i/3(u) 


7re 0 c 2 a;7(ci; 


-lmG(r, r , oj), 


(42) 


where c is the light velocity and 7 ( 0 ;) has been introduced in appendix B. By choosing 


/3(co) = 2he 0 co 2, y( y co), 


(43) 


we have, 


[A(r,t),-e 0 E(r',t)\ = — [ d(n^lmG(r, r, u). (44) 

tt J o c 2 

By considering Eq.(62) and Eq.(63) in appendix A, and using the polar coordinate c 0 = 
Itale*^, the right hand side of above equation is written as follows, 



CO 


den—ImG(r, r ,co) = 



2 limi^i-Kx, 


co 


dip—G(r,r ,u) 


n Xf 

-5(r-r 


(45) 


By substituting (45) in (44) the explicit form of canonical commutation relation which is 
identical with Eq.(20), would be obtained, 


[. A{r , t ), —e 0 E(r ,t)] = ihS(r — r). 


(46) 
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C. The fluctuation relations 


According to the relation (16) and the vector potential operator in (36), the electric field 
can be written, 

E + (r,uj) = iu:A + (r,ou) = /j, 0 C'x 

{(ix + + (ix — ^^-z)e ( '~ UoZ ^Q(z)}x 

Vi v 0 

[cl r (x,uj) + cl l (x,uj)}. (47) 

Regarding to the relation (22) and appendix B, the electric field fluctuation can be 
achieved, 

(0| E(r,u)E(r ,u )|0) = 2hu) 2 nolmG(r,r ,cj)S(u — u ). (48) 

On the other hand from (43) one can deduce the form of the a m (u ) and a d (uj ) in terms of 
the media parameters, 

|a m (a;)| = 2An 2 e 0 Ime m (w), 

|a d (cu)| = 2fan 2 e 0 Imed(o>). (49) 

By applying above relations, the noise current fluctuation is given by, 

(0\j(r,u)j(r\u')\0) = S(r - r')S(u - J) 

2huj 2 e 0 (Ime m (uj)<d(—z) + Ime^o;) ©(;?)). (50) 

These results agree with those calculated by the fluctuation-dissipation theorem. 


IV. THE MAGNETIC FIELD VARIATION IN AMPLIFYING AND ATTENUAT¬ 
ING MEDIA FOR COHERENT AND SQUEEZED SPP MODES 


A. Magnetic field 


In order to express the magnetic field by using (16) and (36), the following relations are 


necessary, 


da R {x,uj ) 
dx 

dd R (x, co) 
dx 


ik spp a R (x,U )) + F(x,u), 

-ik spp a L (x,uj ) - F(x,u), 


(51) 
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where 


, , 2k” . i 

FM = fe )5 


dz 


(ix H— —z)e^ ViZ ^ ■ j^ + (x, z,u)Q(—z)0(—z) 


+ (ix 


k 


—z)e { ~ v ° z ) ■ ' d+ 


v o 


j+(x,z,u)Q(z)<d(z 


(52) 


The equations (51) are the well-known quantum Langevin equations. The F(x,lu) is the 
operator of Langevin noise source and satisfies the following commutation relation, 


[F(x, oj), F(x , oj )] —2k 8(x — x )8 (oj — u ). 


(53) 


Therefore, the magnetic field operators of the propagating rightwards SPP modes are, 
H + (r,uj ) = -{^-)^C'a R (x,uj)x 

{{Vi + ^)e^Q(-z) - (z / 0 + ^)e^Q(z)}. (54) 

Vi VO 

On the other hand, the explicit solution for (51) can be derived, 

a R (x,u) = e iks ^ x - x,) a R (x ,u) + J ^ F(y,u)e ikspp{x - y) dy, (55) 

where x > x . This equation connects two operators in different points of space with each 
other. By applying the general property of the noise operator (f(x,cu)) = 0, the average of 
a R (x,uj), at x = 0 , the starting point of rightwards propagation, is given by, 

(a R (x,u)) = e iksppX (a R (u)). (56) 


The average of the megnetic field can be then derived, 

(ff + (r.w)) = -(^)ic'e ,k ^(a R (u))x 

{(*4 + -z) - (vo + (57) 

Vi V 0 


B. Propagation in amplifying and attenuating media 

The ohmic loss property of the metal plays the main role in dissipation of SPP energy 
and reducing the its propagation length. To overcome these problems one uses a gain 
dielectric medium (with doped dye molecules) adjust to the metal medium. By this system 
not only the SPP loss can be compensated but also the SPP modes can be amplified[23-27]. 
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On the other hand in (32), if the imaginary part is positive, then the SPP modes will be 
attenuated. In this condition the amplification in the dielectric can not compensate the loss 
of the metal. In contrast, if it is negative, the SPP modes can be amplified, because the 
gain of the dielectric can overcome the metal loss [28]. In order to investigate the magnetic 
field average (see Eq.(57)), we consider both attenuating and amplifying systems. Moreover, 
because of the quantum nature of (57) we can consider these conditions for two different 
kinds of SPP modes, i.e., coherent and squeezed states [29]. 

1. coherent state of SPP 

In Ref. [29], it is mentioned that the SPP operators, like photon operators, satisfy the 
Bosonic commutation relations. Therefore, we can define the coherent states of SPP as 
eigenvector of its annihilation operator, in accordance with the definition of the coherent 
states of radiation field, 

an(uj)\a) = a\a). (58) 

By considering the attenuating system [28] and the equation (58), the variation of magnetic 
field average (see (57)) in coherent state is illustrated in Fig. (2). 



FIG. 2. The magnetic field propagation of SPP coherent state in an attenuating system. The 
dielectric medium with refractive index nd = 1.3375 — *0.220 and silver with n m = 0.224 + *1.34 at 
A = 350wn are chosen [28]. k" spp = 7.5 x 10' and |a| 2 = 7. 

Fig. (2) shows that for attenuating system, SPP propagation is attenuated because the 
dielectric gain can not compensate the dissipation in metal. On the other hand, for dielectric 
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medium with refractive index rid = 1.3375 — 10.223 the imaginary part of the SPP refractive 
index is negative and k" spp = —2.2376 x 10', so the system can be considered as an amplifying 
system. Fig.3 illustrates SPP magnetic field propagation in the amplifying system. 


„«n - 10 



FIG. 3. The magnetic field propagation of SPP coherent state in an amplifying system. 

Fig. (3) shows that SPP magnetic field along the propagated axis, can be amplified. It 
means that the gain of the dielectric has overcompensated the loss in the metal. 

2. Squeezed states of SPP 

After observing the quadrature squeezing of SPP in a gold waveguide [35], the quantum 
mechanical description of these states have been given [29]. The squeezed state of SPP is 
generated, similar to the squeezed state of radiation field, by applying the squeezed operator 
on vacuum, 

10 = «(0I°) = ex P(^f a 2 - ^a t2 )|0), (59) 

where £ = |£|e j6le , |£| and 9% are the squeezed parameter and compressed angle, respectively. 
If the SPP is prepared in the coherent squeezed state |£,c/), then 

a(u;)|£, a ) = (/ ua'(uj ) — ua'(u))\a ), (60) 

where /j = cosh(|£|) and v = sinh(|£|)e*^. By applying (57) and (60), we can illustrate the 
magnetic field average of SPP squeezed state, as given for an attenuating system in Fig. (4). 

This figure shows that, the propagation of the SPP will undergo more attenuation in any 
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FIG. 4. The magnetic field propagation of SPP squeezed state in an attenuating system. |£| = 1.5, 
0£ = 0, a = \a\e ze and 6 = 1.5 (rad). 


state other than squeezed state, e.g., if the initial SPP state is prepared in squeezed state, 
the amplitude of the magnetic field is greater than the coherent state case (see Fig. 2). 
Likewise, the SPP magnetic field in an amplifying system and squeezed state is shown in 
Fig.5. 


x 10 



FIG. 5. The magnetic field propagation of SPP squeezed state in an amplifying system. 


The comparison between Fig. (5) and Fig. (3) shows that those systems that the SPP modes 
is amplified, the magnitude of magnetic field in squeezed state is several times greater than 
in the coherent state case. 
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3. comparison between SPP coherent and squeezed state in the amplifying system 

In Ref. [36] the importance of the phase on the properties of a squeezed state, has 
been studied. Accordingly we can investigate the influence of the phase by considering the 
variation of the average of SPP magnetic field in (57) and (60) versus 9, where 9 is the 
argument of a! = \a'\e l6 . It is depicted in Fig. (6). 



FIG. 6. The variation of SPP magnetic field average for different phases and different values of 
squeezed parameter, z = 0 and x = 10 nm. 

As is illustrated in Fig. (6), the SPP magnetic field in squeezed state is greater than in 
the case of coherent state. The drastic difference between coherent and squeezed state is 
occurred in 9 = 1.5 Rad. In Fig. (7) we demonstrate the SPP magnetic field for various 
parameters of squeezed states. 



FIG. 7. The variation of SPP magnetic field average for different values of squeezed parameter, 
z = 0 and 9 = l.5Rad. 
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Fig. (7) shows that by choosing an appropriate phase, the magnetic of SPP in squeezed 
state is several times greater than that of coherent state. 


V. CONCLUSION 

The quantization of SPP by a new method, based on the Green’s tensor, provides the 
possibility to distinguish the behavior of different modes of the quantized field, such as co¬ 
herent and squeezed states of SPP. The propagation of different modes of quantized SPP can 
also be investigated in amplifying and attenuating media. We have seen that the behavior 
of these two modes in the said media are distinguished drastically. 

We think, by extension of this method, the quantization of SPP in different geometries, 
other than the plane interfaces, can also be done straightforwardly, and accordingly it leads 
to understand the more interesting properties of SPP. 


APPENDIX A 

The Green’s tensor satisfy the maxwell’s equation when the source term is replaced by 
Dirac delta function( see Eq. (19)): 

, L0 2 

— V x V x G(r, r ,ui) H—-e(r, uj)G(r, r , oj) = IS(r — r). (61) 

c 2 

For large frequency, e(r, co) —> 1 and from above equation one can deduce that: 

co 2 

lim —G(r,r ,co) = —<5(r — r ) (62) 

|w|—■»oo C 2 

The Green’s tensor has also general properties such as: 

G T (r,r ,co) = G(r ,r,u) 

G*(r, r',uj) — G(r, r, — u) (63) 

On the other hand for a nonmagnetic media,by some algebra, one can obtain a very usefull 
relation [37] 

r , c 2 , 

/ dslme(s, u)G(r, s, u) ■ G*(s : r ,u) = —Im G(r,r ,u) (64) 

J UJ 2 
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APPENDIX B 


The right hand side of Eq.(41) can be express versus Green’s tensor. By considering 
Eq.(23) and the Green’s tensor for a metal-dielectric geometry Eq. (33), one can prove that: 


dslme(s, oj)G(r , s, oj) ■ G* ( s , r , to) = 


2 pik S pp\x x | p \x x | 

\ 2{ — k — + —* — }x 
A r\j fVeTiT) fly 


I*. 


spp • v spp 

[(ix + ^^-z)e^ iZ ^Q{—z) + (hr — -^-z)e ( '~ U0Z ' l Q(z)\ x 
G G) 

[(—ix + -^-z)e v t z Q(—z) + (— ix —-^;$)e _l '° 2 3 4 5 6 7 @(z)], 




ZY 


(65) 


where 


tM = (1 + 


I kspp | \ Ime m , | k S pp | N Imerf 


+ (1 + 


■): 


K| 2 Vi + V? M 2 'z'o + z'o 

Equating Eq.(41) and Eq.(65) and considering Eq.(64) leads to: 

,2 of, ,\ Z2 


[A(r,t),~e 0 E(r',t)\= j dw^^^^hnG(r, r,o;). 

J 7T7(a;) OR 
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